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Abstract 

We consider a classical hydrogen atom in a linearly polarized electric field of slow changing 
frequency. When the system passes through a resonance between the driving frequency and the 
Keplerian frequency of the electron's motion, a capture into the resonance can occur. We study 
this phenomenon in the case of 2:1 resonance and show that the capture results in growth of 
the eccentricity of the electron's orbit. The capture probability for various initial values of the 
eccentricity is defined and calculated. 
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I. INTRODUCTION 



Dynamics of highly excited (Rydberg) atoms in microwave fields has been a subject of 
extensive research during the last thirty years. After experiments of Bayfield and Koch 

Qn 
and theoretical work of Leopold and Percival [2j, it was realised that certain essential 

properties of the dynamics of Rydberg hydrogen atoms can be described in the frames of 

classical approach. 

One of the classical ideas in this area is to control the Keplerian motion of the electron 
using its resonant interaction with a wave of slowly changing frequency. Dynamical problems 
of this kind were studied in for a 1-D model and [4] for a 3-D model. In particular, in 
the latter work an hydrogen atom in a linearly polarized electric field of slowly decreasing 
frequency was considered. It was shown that at a passage through 2:1 resonance (i.e., when 
the driving frequency is twice as large as the Keplerian frequency) the system with initially 
zero eccentricity of the electron's orbit is captured into the resonance. In the captured 
state, the electron's Keplerian frequency varies in such a way that the resonant condition is 
approximately satisfied. In this motion the orbit's eccentricity grows, which may result in 
ionization of the atom. 

In the present work we also consider a 3-D hydrogen atom in a linearly polarized elec- 
trostatic field of slowly changing frequency. We study behaviour of the system near 2:1 
resonance using methods of the theory of resonant phenomena, developed in - (see 
also P)- These methods were previously used in studies of various physical problems includ- 
ing surfatron acceleration of charged particles in magnetic field and electromagnetic wave 



111, 



12| . and classical dynamics in a molecular hydrogen 



[10[, slowly perturbed biUiards 
ion In the present paper, we show that capture into the resonance necessarily occurs 
not only in the case of zero initial eccentricity but also if the initial eccentricity is not zero 
but small enough. Moreover, at larger values of the initial eccentricity the capture is also 
possible. Following the general approach, capture into the resonance in this case can be 
considered as a probabilistic phenomenon. We define and evaluate its probability. The ob- 
tained results can be used to broaden the applicability of the resonant control methods for 
Rydberg atoms. 

The paper is organized as follows. In Section 2, we use standard techniques of classical 
celestial mechanics and the theory of resonant phenomena to reduce the equations of motion 
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near the resonance to the standard form. We consider two different cases: one of small 
eccentricity and the other of eccentricity of order 1. In Section 3, we study the case of small 
eccentricity. We apply relevant results of and find the region of so called "automatic 
capture" into the resonance at small eccentricities and calculate probability of capture at 
larger values of initial eccentricity. Section 4 is devoted to the capture phenomenon at values 
of eccentricity of order 1. We calculate the capture probability in this case too. In the both 
cases, the capture significantly changes the eccentricity of the electron's orbit and may lead 
to ionization. In Conclusions, we summarize the results. 



II. EQUATIONS OF MOTION NEAR THE 2:1 RESONANCE 

We study dynamics of a classical electron in a hydrogen atom perturbed by a harmonically 
oscillating electric field of small amplitude /i, linearly polarized along the Z-axis. This system 
is described with Hamiltonian 

n = no + Zficos"^. (1) 

Here Tio is the unperturbed Hamiltonian of motion in the Coulomb field and is the 
perturbation phase. Introduce the perturbation frequency u = /dt. We assume that 
uj = uj{et), < £ <^ 1, i.e. that this frequency slowly changes with time. For brevity, we 
put the electron mass and charge to 1, and use dimensionless variables. 

The unperturbed trajectory of the electron is an ellipse with eccentricity e, semimajor 
axis a, and inclination i. It is a well-known fact from celestial mechanics that the so-called 
Delaunay elements L, G, H, I, g, h provide a set of canonical variables for the system under 
consideration (see, e.g. |14;]). The Delaunay elements can be defined as 

L = v^, G= v/a(l -e2), H = ^/ a{l - e^) cos i; (2) 

/ is the mean anomaly, g is the argument of the periapsis, and h is the longitude of the 
ascending node (see Figure HJ. 
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FIG. 1: The Keplerian ellipse is shown with the bold line. The rest system of coordinates is 
(XYZ), the z-axis is orthogonal to the plane of the orbit. The periapsis and the ascending node 
are denoted with P and N accordingly. 

In these variables, Hamiltonian (Q) takes the form (see ^): 

1 



n 
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Hq + fiHi , Ho = 



sini 

k=l 
bk — Clk 



'sm{kl + 5( - ^) + sm(A;/ + g + '^)) + 



(3) 



;sin(A;/ - 51 - \1>) + sm{kl - 51 + 



where 



'^Xike), bk 



2a 

ke 



Jk{ke) 



(4) 



Here Jki,-) is the Bessel function of integer order fc, and is its derivative. 

In order to avoid possible singularities at e = 0, we make a canonical transformation of 
variables (L, G, H; I, g, h) ^ (P3, P2, Pi, Q3, Q2, Qi) defined with generating function P^i^l + 
g + h) + P2{g + h) + P^h. The new canonical variables [called Poincare elements of the first 
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kind] are expressed in terms of the old ones as follows: 



Pz = L, Q.^ = l + g + h, 
P2 = G-L, Q2 = g + h, (5) 
P^ = H-G, Qi = h. 
As the perturbation frequency slowly varies with time, the system passes through res- 
onances with the unperturbed Keplerian frequency /. Near a resonance, certain terms in 
expression for Tii are changing very slowly. Consider a passage through the 2:1 reso- 
nance. In this case, after averaging over fast oscillating terms, we obtain the Hamiltonian 
describing the dynamics near the resonance: 

^ = -T742 +/^"(^2,P3) sini sin(2Q3-Q2-Qi-^), (6) 

where we introduced the notation a{P2, P3) = (02 + ^2)/4. 

The resonance is defined by 2Q3 = ^^{'^), t = ^t. It follows from the unperturbed 
Hamiltonian (jS)) that Q3 = 1/P^- Hence, denoting the value of P3 at the resonance as Pr, 
we find: 

Our next step is to introduce the resonant phase. We do this with the canonical trans- 
formation [Pi] Qi) H-i^ ( Jj; ^j), i = 1, 2, 3 defined with the generating function 

W = M2Qs -Q2-Qi-^) + J2Q2 + JiQi. (8) 
For the new canonical variables we have 

J3 = P3/2, 73 = 2Q3 - Q2 - Qi - 

J2 = P2 + P3/2, 72 = Q2, (9) 

Ji=Pi + P3/2, 7i = Qi- 

The Hamiltonian function takes the form: 

H = --r^ + fia{J2, J3) sini sin 73 - uj{t)Jz. (10) 

«J3 

Near the resonant value P3 = Pr we can expand this expression into series. With the 
accuracy of order 0(|P3 — Pj,|^) we obtain the following Hamiltonian: 

^ = ~\ pi '^ + fia{J2, P3) sini sin 73, (11) 



where canonically conjugated variables are P3/2 and 73. Introduce notations P = Psjf = 
73 + 37r/2, J = J2- The Hamiltonian (jllj) does not contain 72, hence J is an integral of the 
problem. Another integral is Ji, corresponding to the fact that Delaunay element H is an 
integral of the original system Q. Coefficient sini in (jllj) should be taken at P = P^- We 
have 

I 1p 

Sinzlp o =\ 1 — - —rr;. (12) 

'^=^'- y (j+p,/2)2 ^ ^ 

From now on, we consider separately two different cases: the case of small initial eccen- 
tricity and the case when initial eccentricity is a value of order one. Let us start with the 
first case. 

Assume the initial value of eccentricity is small, though not necessary zero. From Q, we 
have 

Small eccentricity implies that P/2 — J ^ 1. As the system evolves near the resonance, 
small variations of P are essential when calculating the term (P/2 — J) in (jl3|) and less 
important in the other terms. Hence, in these latter terms, we can put P = P^. We write 

J+^P= (^^- J + 2J + P^ ^2(^7 + :^ 1 ^2P, 

where we have used that P/2 ^ J. Thus, we obtain 



aUP)^y^py^^l^J (14) 



and the following expression for the Hamiltonian: 



Introduce so-called Poincare elements of the second kind: 

X = VP — 2Jcosv5, y = VP — 2J simp. (16) 
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The transformation (P/2, (p) 1— > (x, y) is canonical with generating function Wi = ^ cot ip — 
Jif. Change the sign of J^i and, to preserve the canonical form, the sign of x. Thus, we 
obtain the Hamiltonian in the form: 

^1 = I ^""^ - 3^^(P. - 2J) + f^AiH, J, Pr)x, (17) 



where 



p3/2 / rT2 



Note, that if the eccentricity is small, then — 2J 1, and hence the slow variation of 
Pr is essential only in the second term in ()17|1 . In the other terms, Pr can be assumed to 
be constant, say, Pr = Pro, where P^o is the value of Pr aX t = 0. Now, we renormalise the 
Hamiltonian J^i ^P^^J^i to transform it to the following standard form studied in Qj: 

jr^ = (x^ + yY ~ X{x^ + y^) + fix (18) 

with A = 2{Pr — 2 J) and fl = fiA{H, J, Pro). We describe the dynamics defined by Hamilto- 
nian (fTHj) in Section 3. 

Return now to equation (fTT|) and consider the case when the eccentricity is not small: 
e ~ 1. In this case we can put P = Pr when calculating a. Thus we obtain 

^ ^ ^2 = ~f + fiB{H, J, Pr) cosyp, (19) 

where 



B{H, J, Pr) = sini a(J, P,) = Jl 



(J + P,/2)2 4 



aVl — 



aj^{2e) + — J2{2e) 



(20) 



and values of a and e are calculated at P = P^: 



e.^l-(^. (21) 

The system with Hamiltonian function is a pendulum with slowly changing parameters. 
We study the dynamics in this system in Section 4. 



III. CAPTURE INTO THE RESONANCE AT SMALL VALUES OF ECCEN- 
TRICITY 

Dynamics in the system with Hamiltonian function (jlSp was studied in details in 1^. In 
this section we put forward the results of relevant to our study. 

In ffT8|) /t > is a constant parameter, and A is a slow function of time, A ~ e. Assume 
that A > 0. 
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FIG. 2: Phase portraits of the system at fixed values of A; /i = 0.01. a) A = 0.05 < A*, b) 
A = 0.069624 ^ A*. 

On the phase plane {x,y), the values \/P — 2J and (p (see (HSI)) are polar coordinates. 
Note that eccentricity in the original problem is proportional to a/P — 2 J. 

Parameter A is changing slowly, and as the first step we consider the problem at fixed 
values of A. Phase portraits at different values of A are presented in Figures |2l El At 
A < A* = (Figure IlK) there is one elliptic stable point A on the portrait. At A > A* 

(Figure 13 there are two elliptic stable points A and B, and one saddle point C. Separatrices 
/i,/2 divide the phase plane into three regions Gi,Gi2,G2- In Figure Eb, the portrait at 
A = A,, is shown. 

The coordinates of point C are x = Xc,y = 0, where Xc = Xc{X) is the largest root of 
equation 

— \' ' ' = 4x^ - 2Ax + u = 0. (22) 

ox 

At A > A=K, introduce Fc = Fc{\) = J-'i{xc, 0, A) and jF(x, y, A) = jF(x, y, A) — -Fc(A). In G12 
we have JF < 0, in Gi and G2 ^ > 0, on the separatrices JF = 0. 

As parameter A slowly grows with time, A ~ e, curves hyh defined with JF = slowly 
move on the phase plane. On time intervals of order their position on the phase plane 
essentially changes, together with the areas of G12, 6*2- On the other hand, area surrounded 
by a closed phase trajectory at a frozen value of A is an approximate integral [adiabatic 
invariant] of the system with slowly varying parameter A. Therefore, a phase point can cross 
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FIG. 3: Phase portrait of the system at A = 0.1 > A,,,; fl = 0.01. 

li, I2 leaving one of the regions Gi and entering another region. 

Denote with {x(t),y{t)) a phase point moving according to (fTS|) . Without loss of gen- 
erality assume that A = A* at t = 0. The initial point (x(0),?/(0)) can be either inside 
^1 [-^0 = •^(^^(O), ?/(0), A*) < 0] or out of li [jFg > 0]. The following assertion is vahd. 
All points lying inside /i at A = A* except, maybe, of those belonging to a narrow strip 
—kiE^^^ < ^0 < 0, where ki is a positive constant, stay in G12 at least during time intervals 
of order e^^. 

This result is due to the fact that the area of G12 monotonously grows with time, and 
conservation of the adiabatic invariant makes a phase point go deeper and deeper into this 
region. A point captured in G12 rotates around point A= {xa,0), where xa is the smallest 
root of equation ()22|1 . As time grows, \xa\ also grows and point A on the portrait slowly 
moves along x-axis in the negative direction. Therefore, the motion is a composition of 
two components: fast rotation along a banana-shaped curve surrounding A and slow drift 
along X-axis. The area surrounded by each banana-shaped turn is approximately the same 
and equals the area 5(A*) surrounded by the trajectory passing through (x(0),?/(0)) at 
A = A=K. Hence, the average distance between the phase point and the origin slowly grows, 
corresponding to the eccentricity growth in the original problem. 

In yi , it was shown that a point having zero initial eccentricity necessarily undergoes the 
eccentricity growth. The formulated result implies that this is also valid for all the points 
initially [i.e., at A = A*] inside li, except, maybe for a narrow strip close to /i. A typical 
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linear size of this domain is of order /i^/^. In [l5| this phenomenon was described and called 
"automatic entry into libration". 

Consider now the case when the point {x{0),y{0)) is outside li: JFq > 0. With time the 
area inside li grows, and at a certain moment the phase trajectory crosses li. In the adiabatic 
approximation the area surrounded by the phase trajectory is constant: 5(A) = 5(A*). 
Hence, in this approximation the time moment of crossing li can be found from equation 
iSi(A) = iS(A*). Here iSi(A) is the area inside li as a function of A and A is a value of the 
parameter at this moment. After crossing, there are two possibilities: (i) the phase point 
can continue its motion in Gu during a time interval of order at least e^^ [this corresponds 
to capture into the resonance and growth of the eccentricity] ; (ii) without making a full turn 
inside G12 the phase point can cross I2 and continue its motion in G2 [this corresponds to 
passage through the resonance without capture]. The area of G2 also monotonously grows 
with time, hence such a point can not cross I2 once more and return to G12 

It is shown in ^ that the scenario of motion after crossing li strongly depends on ini- 
tial conditions (x(0), ?/(0)): a small, of order e, variation of initial conditions can result in 
qualitatively different evolution. If the initial conditions are defined with a final accuracy 
(5, e <^ 5 ^ 1, it is impossible to predict a priori the scenario of evolution. Therefore, 
it is reasonable to consider capture into G12 or G2 as random events and introduce their 
probabilities. 

Following Q], consider a circle of small radius S with the centre at the initial point Mq. 
Then the probability of capture into Gu is defined as 

V = lim lim (23) 

sk 

where S^^^ is the measure of the circle of radius 6 and S12 is the measure of the points inside 
this circle that are captured finally into Gu- 

Let A = A be the parameter value at the moment of crossing li in the adiabatic approxi- 
mation. The following formula for probability P is valid: 

P=^^^, where/i(A) = -jf ^dt, /2(A) = -£^dt, (24) 
and the integrals /i, I2 are calculated at A = A. Calculating the integrals, one finds P] 

Ji(A) = ^(27r-e), /2 = |, e = arccos - 2^ (25) 
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Here is the angle between the tangencies to /i at C, < O < it. 

Geometrically, formula can be interpreted as follows. In a Hamiltonian system, phase 
volume is invariant. As parameter A changes by AA, a phase volume AV12 enters the region 
Gi2- At the same time, a volume AV2 leaves this region and enters ^2- The relative measure 
of points captured in G12 is {AV12 — AV2)/AV2. The integral Ji in is the flow of the 
phase volume across li, and I2 is the flow across 12- Therefore, V gives the relative measure 
of points captured into G12. 

Note that, rigorously speaking, there also exists a set of initial conditions that should 
be excluded from consideration. Phase trajectories with initial conditions in this set pass 
very close to saddle point C, and the asymptotic results of cannot be applied to them. 
However, this exclusive set is small: it is shown in 1^ that its relative measure is a value of 
order ^/e. 

In Figure E^, capture into the resonance is shown. First, the phase point rotates around 
the origin in region Gi, then it crosses h, enters region G12 and continues its motion in this 
region. In the course of this motion, the average distance from the origin grows, correspond- 
ing to the growth of the eccentricity. In Figure EJd, all the parameter values are the same 
as in Figure 4a, but initial conditions are slightly different. In this case, after crossing li, 
the phase point crosses I2 and gets into region G2. This is a passage through the resonance 
without capture. 

Summarizing the results of this section, we can say the following, (i) Capture into the 
resonance in the considered case results in growth of the eccentricity of the electron's orbit, 
(ii) On the phase plane around the origin [e = 0], there exists a region of size of order fl^^^ 
such that all phase trajectories with initial conditions [i.e. at A = A*] in this region undergo 
a capture into the resonance with necessity ["automatic capture"], (iii) If initial eccentricity 
is larger, and the initial point on the phase plane is out of the region mentioned above, there 
is a finite probability that the phase trajectory will be captured into the resonance. This 
probability is given by 

IV. CAPTURE INTO THE RESONANCE AT ECCENTRICITY OF ORDER 1. 

If initial eccentricity is a value of order one, dynamics in a y^-neighbourhood of the reso- 
nance is described by Hamiltonian (|19p. In this Hamiltonian, = Pri^) is a monotonously 
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^ I* 

FIG. 4: Passage through the resonance on (x, ?/)-plane and corresponding variation of eccentricity 
(bottom plots, Ecc denotes (x^ + y^)^^^); jl = 0.01, £ = 0.01. a) Passage with capture into the 
resonance; b) passage without capture; initial conditions in the cases a) and b) are different. 

increasing function of the slow time r = st. At a frozen value of Pr, this is a Hamilto- 
nian of a pendulum, with elliptic points aX P = Pr, f = 0mod27r and hyperbolic points at 
P = Pr, (f — 7rmod27r. Denote the value of the Hamiltonian at hyperbolic points with Fg. 
The separatrices connecting the hyperbolic points with (p — — tt and (f — tt are defined by 
equation — Fs or: 

P-Pr^ ±^/2/3P^ ^/lJiB{H,J,Pr) Vl + cos(/?. (26) 

The separatrices divide the phase cylinder P G M, G into domains of direct rotation, 
oscillation, and reverse rotation. The area S of the oscillatory domain is proportional to 
Pr\/JiB; it can be shown that it is a monotonically growing function of Pr- 

Now take into consideration the slow growth of Pr with time. This growth produces 
slow motion of the phase portrait upwards on the phase cylinder. At the same time the 
area between the separatrices slowly grows. As a result, phase point initially above the 
upper separatrix, i.e. in the domain of direct rotation, can cross the separatrix and be 
either captured into the oscillatory domain (this is a capture into the resonance), or pass 
through to the domain of reverse rotation. The final mode of motion strongly depends on 
initial conditions: a small, of order e, variation in them can result in qualitatively different 
final mode of motion. Thus, like in the situation described in Section 3, capture into the 
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oscillatory domain in this problem can be considered as a probabilistic phenomenon. The 
probability of capture can be found as follows. As parameter Pj. changes by a small value 
APr, phase volume AV = {-^^ + L)APr, where L is the length of the separatrix, crosses 
the upper separatrix. In the latter expression, the first term is due to the growth of the 
area of the oscillatory domain, and the second term is due to the slow motion of the upper 
separatrix on the phase portrait. At the same time, phase volume AV^ = ^-AP^ enters the 



oscillatory domain and stays inside of it. Hence, the probability of capture can be evaluated 

as 

Straightforward calculations give: 

V = ^ ' =, (28) 



^p,m + ^i + fiP^K 

where K and K are bounded functions of P^ (we do not write the full expressions here 
explicitly for brevity). It can be shown (see |s|) that asymptotically as e (j2ZI) gives 
the correct value of V defined as in with 5*12 denoting the measure of captured points 
inside a small circle of initial conditions. At fi <^ Pr^-, we find from ()28|) that V ~ ^/Ji. 

For a captured phase point value of P remains at a distance of order ^JJi from P^ on 
time intervals of order e~^. Therefore, as it follows from expressions ()21|). as Pr grows, the 
eccentricity along the captured phase trajectory tends to v^/2, and the semimajor axis of 
the electron's orbit tends to infinity. At J > 0, the eccentricity is always smaller than \/3/2, 
and at J < it is always larger than -\/3/2. Note, however, that in the original system, 
as it follows from ((7j), the rate of variation of Pr is large at large values of Pr- Therefore, 
strictly speaking, the asymptotic methods used in this section are not applicable in the limit 

Pr CO. 

In Figure EK, capture into the resonance is shown. In the beginning, the phase trajectory 
encircles the phase cylinder at approximately constant initial value of P. Then the phase 
point crosses the upper separatrix of the pendulum (fT^ . and enters the oscillatory domain. 
Since this moment, the phase trajectory does not encircle the phase cylinder and the average 
value of P grows. The eccentricity also grows. Figure shows a passage through the 
resonance without capture. The phase point does not stay in the oscillatory domain, but 
crosses the bottom separatrix and enters the domain of reverse rotation and continues its 
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FIG. 5: Passage through the resonance on the (P, ip) phase cy Under and corresponding variation of 
the eccentricity (bottom plots, Ecc denotes 1 — — p-{'^'^ ) ; fi = 0.1, £ = 0.001, J = 0.1, = 0.2. 
Initial conditions are different in the cases a) and b). 

motion at approximately constant new value of P. In this case, the eccentricity undergoes 
only a small variation. 



V. CONCLUSIONS 



We have shown that if the frequency of the driving field slowly decreases, there always ex- 
ists a certain probability of capture into the resonance. A capture results in strong variation 
of the electron orbit's eccentricity, and may lead to ionisation of the atom. The resonant 
capture mechanism is a good tool for control of behaviour of Rydberg atoms. Note, that even 
if the capture probability is small (as in the case considered in Section 4), the phenomenon 
is still important. Consider, for example, an ensemble of Rydberg atoms with various initial 
eccentricities in the case when the driving frequency changes slowly periodically. Then, after 
large enough number of these periods, a relative number of order one of the atoms undergo 
the capture. If the capture probability is a value of order ^///, it will happen after 0(/x~^/^) 
periods, which needs time of order e^^/i^^^^. 
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